Small-q phonon mediated singlet and chiral spin triplet superconductivity in LiFeAs 
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We report fully momentum dependent, self-consistent calculations of the gap symmetry, Fermi 
surface (FS) anisotropy and T c of superconducting (SC) LiFeAs using the experimental band struc- 
ture and a realistic small-q electron phonon interaction within the framework of Migdal-Eliashberg 
theory. In the stoichiometric regime, we find the exact s++ gap as reported by ARPES. For slight 
deviations from stoichiometry towards electron doping, we find that a chiral triplet p x + ip y state 
stabilizes near T c and that at lower temperatures a transition from the triplet to singlet s± SC takes 
place. Further doping stabilizes the chiral p-wave SC down to T=0. Precisely the same behavior 
was observed recently by NMR. Our results provide a natural and universal understanding of the 
conflicting experimental observations in LiFeAs. 

PACS numbers: 74.20.-z, 74.20.Rp, 74.70.Xa 



In iron arsenide superconductors are observed not only 
the higher critical temperatures after those observed in 
cuprates but also a plethora of exciting phenomena yet to 
be understood. Perhaps the most challenging compound 
is the stoichiometric LiFeAs [l| which exhibits a supris- 
ingly exotic superconducting phenomenology. Measure- 
ments of NMR 0], inelastic neutron scattering Q and 
specific heat Q in this compound have been interpreted 
in terms of multigap singlet s± SC [EJ. However, the gap 
seen by ARPES indicates an anisotropic, singlet and one 
sign s++ order parameter On the other hand, there 

is also consistent evidence for p-wave spin triplet SC in 
LiFeAs, first reported by quasiparticle interference (QPI) 
[H and NMR [§]. Subsequent high magnetic field mea- 
surements have found that this state can also be induced 
by the field and it is chiral fiol ] . 

The NMR data have posed yet another extraordinary 
puzzle. Depending on small changes of the stoichiom- 
etry, some of the samples exhibit singlet SC, while the 
slightly more electron doped samples show triplet SC @ . 
Moreover, it was recently reported that some samples ex- 
hibit a temperature induced singlet to triplet transition as 
temperature rises [nj . Clearly, not only do we observe in 
LiFeAs the highest temperature, by an order of magni- 
tude, triplet SC state (T c =16-18 K), but in addition this 
state emerges away from any magnetically ordered phase 
and exhibits some of the most surprising phase transition 
phenomena in a SC regime. Therefore, understanding 
theoretically the complex and seemingly conflicting SC 
phenomenology of LiFeAs is arguably one of the most 
challenging issues in the field of superconductivity. 

The electronic structure of LiFeAs is characterized 
by reduced nesting properties and a Van Ho ve p oint 
(VHp) at the center of the Brillouin Zone (BZ) [12J (see 
Fig. 1(a)). Although antifcrromagnctic spin fluctuations 
(SFs) are generally weak in this material [H, HSU) their 
importance for SC is yet unresolved jljj ]. On the con- 
trary, there is clear evidence of an enhanced electron- 
phonon interaction (EPI), strong enough to account for 
lifetime effects and the order of magnitude of the mea- 



sured T c [15j. So far, it has been ar gued that ferromag- 
netic SFs may drive the triplet SC [16|. On the other 
hand, the s++ SC has been interpreted in terms of or- 
bital fluctuations assisted by phonons [17j . Even if either 
mechanism is proven relevant for LiFeAs, it still remains 
unclear how it could explain coherently the conflicting 
gap symmetry observations. 

In this Letter we demonstrate that all of the aforemen- 
tioned exotic phenomenology is consistently explained, 
in a unified way, if the enhanced EPI in this material 
has a momentum dependence peaked at small-q, simi- 
lar to the one observed recently in cuprates [L8|. Our 
analysis is based on fully anisotropic, self-consistent cal- 
culations of the SC gap within the Eliashberg approach, 
assuming a realistic small-q EPI and using the ARPES 
resolved multi-band structure for LiFeAs as reported by 
the experiments. Remarkably, using an overall electron- 
phonon coupling (EPC) strength of A = 1.38 consistent 
with phonon damping effects in the ARPES spectra and 
reasonable characteristic phonon frequencies of about 100 
K, the same Eliashberg calculations produce the right 
critical temperatures as well. 

Small-q EPI can occur due to strong ionic and/or 
Coulombic effects l3"27|. Note that large ionic polar- 
izabilities that may lead to enhanced dielectricity have 
indeed been reported in pnictides 28, 29|. It was early 
on understood that small-q EPI may lead to unconven- 
tional SC of d-wave type in the cuprates [25[ and recent 
findings support the relevance of this type of interac- 
tion in these materials [3, 27 1. Other unconventional 
SC materials, including k-BEDT organic salts [3(|, the 
heavy fermion UPd2Al3 [3l[ and cobaltites [32| have also 
been shown within a BCS framework to be compatible 
with this generic picture. Recently, the first multiband 
BCS approach with the small-q EPI mechanism in pnic- 
tides [33(| established that the reported unconventional 
SC states are compatible with the measures of finite iso- 
tope effect [HI . Our present fully momentum dependent 
Eliashberg calculations specifically dedicated to LiFeAs 
are of unprecedented precision. The quality with which 



2 



E7" 

a) 



r\ m a ' 



b) 



/i— lOmcV : 



FIG. 1. (Color online) (a)-(c) The FS of LiFeAs colored by ve- 
locity for /u=[0, 10, 50]meV. The a pocket is very shallow and 
posseses high and isotropic DOS. The remaining FS pockets, 
exhibit DOS maxima along F-M that are enhanced by doping. 



the experimental behavior of LiFeAs is reproduced as well 
as the fact that a chiral triplet SC state in a real material 
is plausibly associated with a phononic mechanism, has 
broad implications for the whole field of unconventional 
superconductivity. 

We parametrize the EPC as an effective, yet realis- 
tic, kernel separable over momentum and Matsubara fre- 
quency: A(k, k'; n, n')=\\ s j c i \ nn i , where X nn > has an Ein- 
stein form and Ak.k' = g^k^k'l 2 ' w ith V p h an effective 
potential and No the DOS on the FS. The momentum 
cutoff q c selects the small wave vectors in the attrac- 



tive phonon part while at larger wave vectors the repul- 
sive Coulomb pseudopotential may prevail. For example, 
such a cutoff may arise naturally due to the interplay be- 
tween enhanced dielectricity and the long range nature 
of the Coulomb potential. The simplest approximation 
to this case maps q c to the Thomas-Fermi (TF) screen- 
ing wavevectqr, Qtf oc \J Nq/e where e is the dielectric 
constant 18j, |2CJ, [21[. In the case of LiFeAs, this pro- 
vides with the rough estimate that q c should decrease 
away from stoichiometry, i.e. from the VHp. A small 
q c leads to a momentum decoupling (MD) situation [2~ij ]. 
In the MD regime the gap function loses its rigidity in 
momentum space and becomes correlated with the vari- 
ations of the DOS, leading to FS momentum dependent 
and possibly unconventional SC. For a multiband system, 
this also leads to the dominance of the EPI for intraband 
channels, while the Coulomb pseudopotential provides a 
repulsive interband coupling, allowing for a sign alternat- 
ing SC gap [33|. Since by definition: A = (A(k, k'; 0))££„ 
in the following we fix A to the experimentally reported 
value [HJ by adjusting V p h, so that for each value of q Cl 
A = (A k ,k')k,k' = L38 > where (. . .)£ s is a FS average. 

For the band structure of LiFeAs we use a four band 
Tight-Binding (TB) model elaborated by the IFW group 
as a fit to ARPES data [HI, ISII . Hence our realistic in- 
put, captures the experimentally resolved FS of LiFeAs 
extremely well (FigOJa)). We model the effect of doping 
in the rigid band approximation by substracting a chem- 
ical potential, /x, from all four bands. For fi > (// < 0) 
the system is electron (hole) doped. We refer to the inner 
hole, outer hole, inner electron and outer electron pocket 
as a,/3,7 and 5, respectively. In FigQJa)-(c) the FS of 
this model is shown, colored by velocity at different val- 



ues of n . The presence of a VHp at T is evident by the 
blue color denoting a very high and isotropic DOS over 
a. This pocket is very shallow thus, a slight deviation 
from stoichiometry, which translates to /i > 5.2meV, is 
enough to remove it from the FS as shown in Fi^IJb)-(c). 
The remaining pockets exhibit some anisotropy, which is 
enhanced with doping, and maximum DOS along T-M. 

Having analyzed the input of our theory, we now dis- 
cuss briefly (for details see the SOM) the technique that 
we employ to obtain our results. Our starting point is the 
fully anisotropic Eliashberg equations (e.g. see ref.(3(J). 
The EPC that we use implies that the SC gap function 
is also separable and can be written as: n = A. n g^, 
where A„ is the gap amplitude distributed over the Mat- 
subara frequencies and g^ contains the momentum de- 
pendence of the SC gap (|c/k < !)• Utilizing this prop- 
erty, we end up with two eigenvalue equations at T=T C 
instead of one, as in standard Eliashberg theory. The 
first one is the usual eigenvalue problem but extended to 
account for the full anisotropy of the gap: 

N 
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" _ (Igk Pp 3 ' Ak = ( Ak : k ')k' 6 and even frequency pair- 
ing is explicitely assumed. The Coulomb pseudopoten- 
tial, n* , is taken as band independent. The T c can be 
calculated by Eq.fll]) provided that g k is known. The lat- 
ter can be obtained by the second eigenvalue equation: 



^((Ak.k' -M*)ffk'} 
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k' 



A9k 



(2) 



with Zk = 1 + A k . Eq. ([2]) , when solved self- consistently , 
yields the symmetry and the exact momentum depen- 
dence of the SC gap. This modelling, permits us to cal- 
culate accurately <?k with a large resolution of grid points 
and then obtain the exact T c within the full Eliashberg 
framework. The linearized equations do not provide in- 
formation on the SC gap below T c , e.g. whether the SC 
symmetry changes with T. In order to estimate such a 
possibility, we also calculate g^ at T=0 using the BCS- 
like self-consistent equation: 

A k = ^k 1 ((A k , k ' — M*) A k ' sinh- 1 ( WD /|A k <|))£ 5 (3) 

where u>d is the Debye frequency. The above, supple- 
mented with the respective free energy formula, pro- 
vides an approximate description for the gap structure 
at T=0 in the limit where the feedback of the gap on the 
renormalization function is negligible (for details see the 
SOM). Using the solution of Eq.<(S16|> in Eq.®, we can 
calculate the anticipated T c of the T=0 state. A com- 
bination of the T=0 and T=T C results suffices in order 
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FIG. 2. (Color online) (a)-(c) Self-consistently calculated 
phase diagrams as a function of the Coulomb pseudopoten- 
tial li* and the EPI momentum cut-off q c , for A=1.38 and 
Li=[0, 10, 50]meV at T=T C . The y-axis is logarithmic to high- 
light the small-g c region. The background color grade corre- 
sponds to the calculated T c at each (q c ,Li*) value, (d)-(f) Same 
plots at T=0. Here, the color grade is for calculations of the 
anticipated T c from the low-T momentum structure of the SC 
gap, as described in the text. The notation s++, s± and p± is 
for anisotropic one-sign and extended sign-alternating singlet 
s-wave and chiral triplet SC gap symmetries, respectively. It 
is easy to observe that the /i=0 and /i=50meV results for both 
T limits agree very well with each other. Hence, at these dop- 
ing levels, we expect that the symmetry of the SC gap found 
for T = T c is unchanged down to zero T. On the contrary, 
for intermediate dopings (here n=10meV), there is a large 
parameter regime where the low-T gap symmetry is singlet 
s±, while for higher temperatures, we find chiral triplet SC 
(e.g. compare the situation marked by the diamond symbols 
in Fig.(b) and (e)). These findings are in perfect agreement 
with the recently observed singlet to triplet SC temperature 
induced transitions in LiFeAs [Xl|] (see Fig. 4 for H=0). 



to discern possible changes in the SC symmetry with T. 
We solve Eq.@ self-consistently over a 512x512 grid of 
k-points for various initial forms of g^ including every 
allowed symmetry of the D^h group or random values. 
The solution that is kept, maximizes the eigenvalue p. 
Next, the coefficients a,b,c are calculated and plugged 
into Eq.fl]) which is solved for the T c . The cutoff fre- 
quency used is w c = lOwrj, where for wd we use the 



calculated logarithmic frequency oj/„=100K 37|. These 
steps are repeated varying li £ [0, 50mey], q c £ [5x5-, j£\ 
(the lattice constant a=l) and li* £ [0,0.3] while keep- 
ing A=1.38 fixed. The same procedure is followed for 
Eq. ([S16p where the correct solution minimizes the free 
energy. 

The obtained q c -p* phase diagrams (PD) and the cal- 
culated T c for various dopings at T=T C are shown in 
Figj2[a)-(c). In the stoichiometric regime (Figl^a)), sin- 
glet Ai g states dominate the PD. This behavior can be 
attributed to the respective FS topology of LiFeAs shown 
in Figllja). Due to the VHp at T, electrons from the 
isotropic a-band have a pronounced contribution to SC 
pairing thus favoring isotropic gap structures over SC 
symmetries with nodes such as d or p-wavc. For rela- 
tively large q c values, we find a one-sign anisotropic, s ++ , 
SC gap. Decreasing q c or increasing li* the system, in or- 
der to avoid the pair-breaking Coulomb repulsion, enters 
into a sign-alternating nodeless s± state. Even in the 
latter, the symmetry is not pure s±, but a strong s-wave 
component is superimposed. Hence, the T c of both states 
is affected by li* , since b ^ in Eq. (JTJ) . Notice how T c in- 
creases with decreasing q c since then, the EPC becomes 
sharply peaked and tends to become purely intraband, 
maximizing the contribution to SC pairing. 

From the T c dependence, it can be seen that there is 
a reasonable parameter space where our phonemic theory 
produces the experimental T c «18K of LiFeAs. For this 
parameter space (q c > 0.1) we find that in both s++ and 
s± regimes the gap over the a-pockct is the largest and 
isotropic. This is a result of the combined VHp and MD 
effects. The gap on the other FS sheets shows anisotropy 
that varies depending on symmetry. In the s± state, the 
/3-pockct exhibits maxima along the T-X and the electron 
pockets along the X-M direction, respectively. In the s++ 
state, the /3-pocket exhibits maxima along the T-X and 
the electron pockets along the L-M direction. The latter 
features match precicely the FS dependence reported by 
ARPES for this material 0, [?J . For example, for li* = 
0.13 and q c = 0.25 we find s ++ SC with T c wl8K and FS 
momentum dependence as shown in FigO^b). 

In the non-stoichiometric case the a-pocket is removed 
from the FS as can be seen in Figlljb)-(c) for li =10 
and li =50 meV, respectively. The remaining anisotropic 
FS possesses minima along the directions parallel and 
perpendicular to T-X that are enhanced with li. Hence, 
the system tends to allow nodes or gap minima in these 
directions, possibly favoring d xy or p-wave SC. The re- 
spective T=T C PDs are shown in Figl2jb)-(c). Increasing 
li* the s± region shrinks and at q c < 0.1, a p-wave SC 
state (<?k ~ siiik x or smk y ) is indeed stabilized. Since 
in our formalism SC pairing is intraband, this is neces- 
sary spin triplet in order for the antisymmetry of the 
electron wavefunction to be satisfied. Free energy cal- 
culations at T=0, indicate that the stable solution is in 
fact chiral (g^ ~ s'mk x ± isvaky), which we hereafter 
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TABLE I. Summary of key results for A=1.38, ^*=0.13 
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FIG. 3. (Color online) Self-consistently obtained mo- 
mentum dependence of the SC gap (\gk\) at T=T C for 
A=1.38 and u*=0.13 in the (a) stoichiometric s++ regime 
(/i=0meV, g c =0.25) and (b) in the non-stoichiometric p± 
regime (fi—lOmeV , q c =0.032). The gap on each band is nor- 
malized to unity in order to highlight the location of gap min- 
ima/maxima. The s++ SC gap in (a) has the exact FS mo- 
mentum dependence reported by ARPES 0, Qj- 



refer to as p± . The absence of the VHp affects also the 
T c which gradually decreases with \x. For small n, this 
decrease can be compensated with a decrease in g c , as is 
evident in FigJ^b). Remarkably, for (q c ,n*) values that 
yield T c «18K, the system is in the p± state. The oc- 
curence of p± SC only for \x > within our theory is 
in perfect agreement with QPI [i[ and particularly NMR 
0,1m reports, where in the latter it has been explicitely 
associated with non-stoichiomctry. Moreover, it resolves 
the conflict between the reported triplet SC and ARPES 
data that indicate a strong EPI in LiFeAs [l5| . 

We now compare our T=0 results with the ones at 
T=T C . The T=0 PDs for ^=[0, 10, 50]meV are shown 
in FigEtd)-(f). For ^=0 (FigE^a) and (d)), it is easy to 
notice that not only the PDs but also the calculated T c 's, 
in the parameter region of interest discussed above, agree 
nicely. The same also holds for the ^=50mcV results 
(FigHJc) and (f)). In fact, for fi > 50mcV, the two PD's 
become exactly the same. Hence, we conclude that the 
T=0 approximation works very well both in determining 
the PD as well as predicting the T c . On this firm basis 
we estimate that for /i=0 (fj, >50meV) LiFeAs is in a s ++ 
(j>±) state down to zero T. 

Our solutions for < fj, < 50meV are amenable to 
an interesting interpretation. In this regime, we sys- 
tematically find that for relevant (q c ,fj,*) values, the SC 
symmetry at T=0 is s± , while at T=T C it is p± . More- 
over, the expected T c for the s± is smaller than the T c 
of the p± state. Thus, we find clear evidence of a tem- 
perature induced transition from singlet s± to triplet p± 
SC that can take place for small deviations from stoi- 
chiometry in LiFeAs. Remarkably, such transitions have 
been reported very recently by NMR in samples that 
are slightly non-stoichiometric (llj . For example, for 
^=10meV, /U*=0.13 and g c =0.032, we find two succes- 
sive transitions, to p± SC below T Cl «18K and to s± SC 
below T C2 «16K (see diamond symbol in Fig|2jb) and 
(c)), in perfect agreement with the findings of the latter 
experiments summarized in Fig. 4 of ref.[llj. Due to the 
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small-<7 c values in this region, the enhanced MD leads to 
a SC gap structure whose FS momentum dependence fol- 
lows closely the FS DOS. Thus, we predict that the SC 
gap in both the s± and p± regions should exhibit maxima 
along T-M, as is shown e.g. in Figl3{c). Our key results 
for ^i*=0.13 are summarized on Table HI 

Finally, given the close proximity of s++, s± and p± 
states found within our theory, it is natural to expect that 
an applied magnetic field would induce a singlet to triplet 
transition in LiFeAs, either by suppresing singlet pairing 
or by enhancing the FS DOS anisotropy or both, thus 
favoring p± SC. Hence, our mechanism could provide a 
plausible explanation for the experimental reports of a 
field induced triplet SC phase in this material [9I- 11|. 



In conclusion, we have presented a theory from the 
small-q phonon perspective, within which a plethora of 
controversial and yet unexplained experimental findings 
in SC LiFeAs are coherently understood. Based on re- 
alistic, fully anisotropic Eliashberg calculations, we have 
demonstrated how the interplay of this small-q EPI and 
the intrinsic FS properties of LiFeAs, governed by the 
presence of a VHp at the center of the BZ, result in a 
delicate balance between singlet s ++ ,s± and triplet chi- 
ral p x + ip y SC. Furthermore, we explicitly showed that 
slight deviations from stoichiometry and/or changes in 
temperature may favor one of these states. Our results 
resolve the conflict between a strong EPC and one sign 
singlet SC reported by ARPES and unconventional and, 
most interestingly, chiral triplet SC observed by other 
probes. In addition, they provide a systematic under- 
standing for the occurence of triplet SC in non stoichio- 
metric samples as well as the temperature induced sin- 
glet to triplet SC transitions reported very recently. The 
accuracy of the present results and the exotic charac- 
ter of the involved SC states, establish that the small-q 
EPI mechanism should be considered on an equal footing 
with spin and other purely electronic mechanisms, in the 
analysis of any unconventional superconductor. Clearly, 
even the most exotic SC states like the chiral spin triplet 
reported here, should not exclude apriori our phononic 
mechanism. 
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SMALL-Q ELECTRON-PHONON COUPLING AND ITS SEPARABILITY 

In this section, we discuss the simplest situation when a small-q electron-phonon coupling (EPC) may arise and 
show how the obtained EPC maps to a sebarable one without affecting any of the basic physics. The EPC function 
is defined as: 

Ml ,, 1 „,„ 1 = J /-™| 41 (S1) 

j {u) n - uj n ,y + n 1 

where a 2 F(k, k'; 57) is the momentum dependent electron-phonon spectral function (EPSF), which for a single phonon 
branch is: 

a 2 F(k,k';fi) =iVo|<? kk '| 2 8 (fl- (S2) 

with ^kk' the electron-phonon matrix element (EPME) [l[ . Perhaps the most common situation when small-q processes 
dominate can take place when large dielectricity and the long range nature of the Coulomb potential are considered, 
as first discussed by Abrikosov (ref. [20] of main text). As a paradigm, we analyze this case in some detail. Due to 
the modulation of the screened Coulomb potential, there exists a first-order coupling between electrons and phonons 
Q. In this case, the electron-phonon matrix element is given by: 



g « = V 2M^- ■ q] t (S3) 
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where q = k — k' is the exchanged momentum, M is the ionic mass, w q is the renormalized phonon dispersion, e q is the 
phonon polarization, Vq = ^p- is the bare Coulomb potential, e is the static dielectric constant and e q is the static 
dielectric function. The renormalized phonon dispersion is just w q = -^=, where O p is the phonon plasma frequency. 

2 

In the simplest picture, the dielectric function can be written in the Thomas- Fermi approximation as: e q = 1 + ■Of-, 

where qxF is the Thomas-Fermi screening wavevector and qrF = \/4:ne 2 No/e. 
Compining the above relations, the EPC can be written in the following form: 

A(q;m) = N V ep -^- 2 (S4) 

where the effective parameter V ep includes all remaining terms and w rn = Iwi-kT is a bosonic frequency (to = n — n'). 
The EPC has a small-q momentum dependence peaked at q — qtf and a Lorentzian frequency dependence whose 
width is determined again by qxF- Thus, even in this rough approximation, the EPC acquires a small-q momentum 
dependence with the momentum cutoff defined by the Thomas-Fermi screening wavevector. The peak at small-q gets 
more pronounced by reduction of the density of states or an increase in the dielectric constant. At low energies, (w=0 
or T=0), Eq.([S4| reduces to: 

A(q) = N V ep 2 1 2 (S5) 
Q i Qtf 

Increasing the ratio T/O p , Eq. (|S4[) approaches the form of Eq. (|S5|) for m ^ 0. For parameters relevant to LiFcAs this 
is the case near T ps T c (i.e. see FiglSip. Thus, A q m behaves like a separable function over momentum and frequency 
in the physical range that we are interested in. 



Separable small-q EPC 

Let us now assume an effective separable EPC. This can be achieved, for example, if the EPSF is separable and 
the frequency dependent part is described by an Einstein spectrum 0, Q] : 

aiF(k,k';Q) = A kk ^J(fi-0 B ) (S6) 
Inserting Eq. ([S6)l into Eq. ([ST|) one easily gets: 

A(k,k';n,n') = A kk > ? ^2,02 = ^kk'A»„' (S7) 

At low energies (n = n'/T=0) the above reduces to: A kk ' which we assume to have a small-q structure within a cutoff 
q c . Thus, we can write: 



A sep (q;TO) = N V ep ] (S8) 
q -f q c w m -|- Ji £ 

Obviously, the above matches exactly Eq. (|S4[) at low e nerg ies (w—0/T=0) and for q c = qrF, Q P — &e- In our 
calculations for LiFeAs, we set 51 p (£)=cj; og =100 K (ref. 37| of main text). We observe that already at T=10K the 
EPC of Eq. (|ST|) and Eq.JS8| agree very well as it is shown in FiglSTl 



MULTIBAND MOMENTUM DEPENDENT ELIASHBERG FORMALISM 



Our starting point is the system of coupled anisotropic Eliashberg equations at the finite temperature: 



Zw,n — 1 



TTT 



^<A(k,k';n,n') 



A 



k' n' 



\u n i I <u a 



Z^ n A k ^ n = 7t-T ^ ([A(k,k';n,?i') - /i*(o 



\FS 
/k' 



A k , 



IV 



(S9) 
(S10) 
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A(q,5) 






A(q,10) 






— A sep (q,o) : 






- A scp (q,5) 






A scp (q,10) 









t-- . ---^n=-.-.T.v^-.- ? , .-.r.-.^ ra -^^- i -^ 

0.0 0.1 0.2 0.3 0.4 0.5 

q/jT 



FIG. SI. Momentum dependence of A(q, m), A sep (q, m) for g TF ( c )=0.3, f2 p (E) = 100K at T=10K and different Matsubara 
frequencies. The agreement between A and \ aep is evident. 



where the nth Matsubara freuency w n = (2n + 1) 7rT, Zk. n = Z(k, iLo n ) is the strong-coupling rcnormalization param- 
eter, Ak, n = A(k, iui n ) is the SC gap function, A(k, k'; n, n') the EPC and fi*(uj c ) is the Coulomb pscudopotential 
which comes with a cutoff lo c in the Matsubara frequency summation and is taken as band independent. The notation 
(. . .) k ?s = J2 k ffi (. . . ) means an average over the entire FS, where Nq is the Density of States (DOS), = V£k F | _1 

the angularly resolved DOS (arDOS) and £k the energy dispersion at the Fermi level (X^^ ]vo = ^' ^ n tn ^ s 
formalism, all band structure effects are conveniently encoded in the FS averages and the anisotropy of the FS is fully 
taken into account. 

Eq. (|S9[) - (|S10|) provide the strong coupling description of a phonon mediated, singlet or unitary triplet, supercon- 
ductor (e.g. see |5[) having an arbitrary number of bands contributing to the FS [fj. In this sense, Ak,n,Zk,n are 
global quantities defined over the entire FS. For an i-band SC, it is possible to degrade the above into a system of 
2i-coupled equations for A k n , Zt^ defined on each separate band, however by doing so, the momentum dependence 
of the EPC is lost (e.g. see ref.[6[). Thus, throughout this study, we work in the most generalized framework as 
defined above. Note that, since the values of fj, that we consider are very small, we have neglected the energy shift 
self-energy term. 



Equations at T c in the separable model 

Linearizing the above system at T c , introducing the separable approximation: A(k, k'; n, n') = Xk.k'^n,n' and re- 
cxprcssing the Matsubara sum on positive frequencies yields: 

i " 

^k,„ = l + A k — — (A(0) + 2 V A(m)) (Sll) 

m— 1 

UJ n / <CJ C 

4, n Ak,n = ^2 Wn'(|A k ,k' (A„_„' + A n+n / +i ) - 2^*(cj c )|A k ',„')k'' S (S12) 



n'=0 



where we have explicitely assumed even frequency SC and Ak=(Ak.k')k /S ' Wn ~ 2n+i > an< ^ we nave used the relation: 
Yin' ^""'jZTTT = A(0) + 2 Ylm=i A (to). From Eq. (|Sll|) . it is evident that Zv „ has become separable. Inserting Eq. (|Slip 
into Eq. ([S12[ ) we see that the SC gap has also acquired a separable structure and thus, we can write: Ak tn =A n ffk, 
where A„ contains the frequency and <7k| < 1 the momentum dependence, respectively. 

Inserting Eq.JSTTJ into Eq. (|S12[) . multiplying both sides with gt — T , — and then taking the FS average gives: 

\ 1 9k I / k 



A„= £ 



(A n _„- + A„+„,+iJ - Zfi (u c ) , 2)FS J /, |2\FS dnn ' 



(i<? k i 2 )r it ™ ^' (i5 k i 2 >r J (i5ki 2 )^ 



A(0) + 2]T A(m) 



w n >A n < (S13) 
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which, after some rearanging can be written as: 



N 

[a ■ (A(n - n') + A(n + n' + 1)) - 26 - c ■ 5 nn , x (a(0) + 2 ^ A(m)) ' 

n ; — m— 1 



w„ 



A n , = 



with A„/ = A n /«v, a = „ kk ?X'f*' , b = y % 'MJ and c = K k, ,l^ s . The above is an eigenvalue problem 

similar to the one encountered in usual Eliashbcrg theory [7[ but extended to include the full momentum dependence 
of the SC gap and it can be solved following standard methods [f| for the T c provided that <7k is known. Obviously, 
for gk — I the usual isotropic Eliashberg equation is retrieved. 

Since A kj „ = A„<7 k , we can apply a square-well ansatz || to Eq. (|Sll|l - (|S12|l in order to extract an equation for the 
momentum dependent part of the SC gap: 

x/1 ,/ a J Ak,k' for both \u n \, \u n >\ < w c , jA(T)g k , \u n \ < u> c , 

A(k, k ; n, n ) = <^ and A„ 5k = 4 

I (J otherwise II), |w„ > uj c 

It is easy to observe that within this model Z^ n becomes: 

^k = 1 + A k (S14) 

and the gap equation takes the form: 

\w n i I <u c 

Z k5k A(T)= w «'({2A k ^ -2 M > c )}. 9k ,)£ S A(T) 

ra'=0 

which after a standard summation on Matsubara frequencies yields: 



where p —In 



1.13cj c 



^ k 1 ((A k ^-/i*) 5k ,)k' S = P5k (S15) 

rp . This is an eigenvalue equation for the momentum part of the SC gap alone, which when 
solved self-consistently provides the exact g k . Notice that in this case _g k is not just a form factor corresponding to an 
irreducible representation, but it is rather a superposition of form factors and all their harmonics of all the allowed 
irreducible representations, i.e. it is the realistic SC structure favored by the system's specific characteristics. 

Equations at T=0 in the separable model 

At T=0, the Eliashberg equations can be treated approximately to give a closed form, BCS-like equation. This can 
be achieved by neglecting A ki „ in Eq. (|S9|) and further applying the square- well ansatz (see e.g. ref.Q Using 
Zk = 1 + A k in Eq. (|ST0l> and taking the zero T limit (T J2 n , ->• 2 f£° f^) gives: 

A k = Z k 1 ((A k ,k' -M*) Ak'Sinh- 1 f-^JL-W (S16) 



After self-consistently solving the above equation, <7 k is retrieved. The correct solution is determined by minimizing 
the respective free energy difference between the normal and the SC state within Eliashberg theory . A momentum 



dependent expression for the latter is (12J: 



5F = -vrTiVo ^((M + A ln ~ Kl) ( Z ^ ~ Z% n . Kl = )) 



FS 
k 



Applying the square-well approximation, taking the T=0 integral over Matsubara frequencies and for = Z k we 
find the respective condensation energy as: 



SF = N iu c (Z k (Lu c - yJu>* + Al))£ s (S17) 

The isotropic version of Eq. (|S16|) has been previously considered in multiband Eliashberg theories, where it has 
been shown that it, in fact, captures the full Eliashberg results very well (e.g. see ref.[l(J). Our results for LiFcAs, 
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indicate that this approximation scheme works also well, when the full momentum dependence is also included in the 
calculations. 



' varelogi@mail.ntua.gr 

[1] G. Grimvall, The Electron- Phonon Interaction in Metals (North-Holland Publishing Co., Amsterdam), 1981. 

[2] J.R. Schrieffer, Theory of Superconductivity (W.A. Benjamin, New York), 1964. 

[3] M.Daams and J. P. Carbotte, J. Low Temp. Phys. 43, 263 (1981); 

[4] A. J. Millis, S. Sachdev and C. M. Varma, Phys. Rev. B 37, 4975 (1988). 

[5] PB. Allen and B. Mitrovic, Solid State Phys. 37 (1982). 

[6] H.J. Choi, M.L. Cohen and S.G. Louie, Phys. Rev. B 73, 104520 (2006). 

[7] PB. Allen and R.C. Dynes, Phys. Rev. B 12, 905 (1975). 

[8] J.P. Carbotte, Rev. Mod. Phys. 62, 1027 (1990). 

[9] E.J. Nicol and J.P. Carbotte, Phys. Rev. B 71, 054501 (2005). 
[10] O.V. Dolgov et al, Phys. Rev. B 79, 060502(R) (2009). 
[11] J. Bardeen, M. Stephen, Phys. Rev. 136, A1485 (1964). 
[12] H.J. Choi, M.L. Cohen and S.G. Louie, Physica C 385, 66 (2003). 




